v=2.4e-03;
dm=0.0013;
w=10;
a=load('7_bodyforce.txt');
[nx,ny]=size(a);
x=1:nx;
x=x*1000;
y=0.003*(x.*0+1);

D=dm+(w^2)*v^2/(210*dm);
z=D*(x.*0+1);


plot(x,a,'b');
hold on
%plot (x,y,'-r');
%hold on
plot(x,z,'-k');
hold off
